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The gamma-ray excess observed by the Fermi-LAT in the Galactic Center can 
be interpreted by the dark matter annihilation to bb via a light pseudoscalar in the 
NMSSM. It is interesting to note that the corresponding singlet scalar is useful to 
achieve a strongly first order phase transition required by the electroweak baryo- 
genesis. In this paper, we investigate the possibility that the NMSSM model can 
simultaneously accommodate these two issues. The phase transition strength can 
be characterized by the vacua energy gap at zero temperature and be sufficiently 
enhanced by the tree-level effect in the NMSSM. We find that the annihilation of 
Singlino/Higgsino DM particles occurring close to the light pseudoscalar resonance 
is favored by the galactic center excess and the observed DM relic density, and a 
resulting small n/X and a negative A^ can also lead to a successful strongly first 
order electroweak phase transition. 

PACS numbers: 95.35.-hd, 12.60.Jv, ll.lO.Wx 

I. INTRODUCTION 

The weakly interacting massive particles (WIMPs) are the most popular and attractive 
dark matter (DM) candidates. WIMPs with masses of 0{1) — (P(IO^) GeV naturally provide 
the correct DM relic density via the freeze-out mechanism in the early Universe. The gamma 
rays produced by present DM annihilations are very good probe to reveal the microscopic 
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nature and the distribution of DM due to the simple propagation and low energy loss. Such 
DM induced gamma-ray signatures from astrophysical sources with high DM densities, such 
as dwarf galaxies, galaxy clusters and the Galactic Center, have been extensively studied. 

Recently, an extended gamma-ray excess of a few GeV in the GC was discovered in the 
Fermi Large Area Telescope data CHm. Although the galactic center excess (GCE) is 
highly statistically signihcant, the diffuse gamma-ray background model in the GC would 
introduce some unclear systematic uncertainties. After considering such uncertainties, some 
studies still conhrmed the existence of the GCE m [TT] . The origin of the GCE may be 
astrophysical sources, such as a population of millisecond pulsars (MSP) |T2[[T3]. However, 
the MSP scenario is not easily explained the energy spectrum and spatial distribution of the 
GCE [m [T^ . Additionally, whether emissions from all MSPs can account for the total flux 
of the GCE is still under debate ra. 

A very attractive interpretation for the GCE is the DM annihilation. The spatial distri¬ 
bution of the GCE is compatible with the square of Navarro-Frenk-White DM distribution 
with a slop ~ 7 = 1.2. The energy spectrum of the GCE is well £t by a ~ 30 — 40 GeV DM 
particle annihilating into bb with a cross section of (av) ~ 2 x 10“^® cm^/s. A ~ 7 — 10 DM 
particle annihilating into rr with a cross section of (av) ~ 5 x 10“^^ cm^/s is also possible 
to £t the GCE, but the quality of this fit is lower than the fit for bb channel [8]. These 
models can be directly tested by the Fermi-LAT observations of dwarf spheroidal galaxies. 
The latest Fermi-LAT results have improved the current limits by a factor of 0{1) and set 
constraints on the DM models accounting for the GCE |T6]. However, considering the un¬ 
certainties from the DM density profiles in the GC and dwarf galaxies, the tensions between 
different Fermi-LAT observations can be relaxed d- 

In literatures, simplihed models containing new mediators connecting DM particles to bb 
have been proposed to explain the GCE [THH2n]. Since the masses of mediators required 
by the GCE are within the region of (P(IO^) — (P(IO^) GeV, they can be directly on-shell 
produced at colliders and then decay into DM particles. At direct detection experiments, 
even there is no directly interaction between the DM particle and light quarks at tree level, 
DM-b quark interactions may also induce DM-nucleon scatterings at loop level. Therefore, 
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the results of collider and direct detection experiments would set constraints on simplified 
models accounting for the GCE [T8H20]. 

In principle, UV complete theories can be mapped into simplihed models [21]. Super- 
symmetry (SUSY) is an very attractive and theoretical motivated new physics model; it can 
provide suitable DM candidates and mediators to explain the GCE. In the SUSY model, 
DM particles can annihilate into bb via t-channel exchange of light sbottoms. The limits 
on directly sbottom pair production from LEP and LHC results are very stringent. A more 
promising s-channel annihilation channel is through a pseudoscalar xx a ^ bh. Compar¬ 
ing with the annihilation mediated by a scalar, this process is not suppressed by the small 
DM velocity in the Galaxy. Since the pseudoscalar mediator accounting for the GCE is 
required to be lighter than ~ 100 GeV, the LHC Higgs searches have set strong limits on 
the pseudoscalar in the context of the minimal supersymmetric standard model (MSSM). 
An economical SUSY extension model evading all collider and direct detection limits is the 
next-to-minimal supersymmetric standard model (NMSSM) [22], which provides a correct 
SM-like Higgs mass ~125 GeV and solves the /i problem in the minimal supersymmetric 
standard model. In the NMSSM, the pseudoscalar mediator is the lightest CP-odd Higgs. 
As pointed in Ref. [2I|, the GCE and DM relic density can be simultaneously interpreted 
by the Singlino/Higgsino DM resonant annihilation or the Bino/Higgsino DM off-resonant 
annihilation via the exchange of light pseudoscalar (see also Refs. [231126] ). 

Another fundamental problem in the Cosmology is the origin of the baryon asymmetry. 
An attractive scenario is the electroweak baryogenesis (EWBG), where a strongly Erst order 
electroweak phase transition (SEOEWPT) is required to avoid the washout of the generated 
baryon asymmetry of Universe (BAU)^ . The EWBG connects the Cosmology to the Higgs 
physics and can be tested at high energy colliders. In the MSSM, a light stop is required 
to successfully trigger the SEOEWPT [2HH3T]. After the Higgs discovery, a big challenge 
is how to simultaneously explain the correct Higgs mass ~ 125 GeV while evading strin¬ 
gent constraints from Higgs measurements. In fact, recent results of LHC Higgs and stop 

^ In fact, SFOEWPT is not enough for the successful EWBG and the strength of the SEOEWPT may affect 
the CP violation source which drives the EWBG ETj. 
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searches have ruled out the mass window of the light stop in the MSSM [321135] . If the MSSM 
Higgs sector (consisting of two Higgs doublets) is enlarged to include an additional singlet 
superheld, we have an alternative supersymmetric framework, namely the next-to-minimal 
supersymmetric Model (NMSSM), which could provide the successful SFOEWPT, as been 
appreciated for a long while |HB1143] . As pointed in Ref. [12], the SFOFWPT would occur in 
Rk ~ —1 and positive < (P(IO) in NMSSM, where = AnVs/Af^ is a critical parameter. 
In this case, the lightest CP-odd Higgs is very light. It is interesting to hnd the connec¬ 
tion between the DM phenomenology and the EWBG. In this paper, we have studied this 
problem in detail and found the suitable parameter space simultaneously accommodating 
the GCE, DM relic density and SFOEWPT in the NMSSM. 

This paper is organized as follows. In Sec. H, we reviewed and resummarized the key 
of the strong Erst order electroweak phase transition in the NMSSM, and we make the 
energy gap analysis which gives the clue of EWPT. In Sec. HI, the interpretations of GGE 
in NMSSM is reviewed, the resonance effect in the situation of k/A 1 is analyzed. In 
Sec.IV, the numerical analysis of the EWPT, GGE and DM relic density are carried out, the 
benchmark scenario which could explain GGE and give rise to the SFOEWPT and correct 
DM relic density is presented. In Sec.V, we summarized the work and give our discussions 
and conclusion. 

II. STRONGLY FIRST ORDER ELECTROWEAK PHASE TRANSITION IN 


THE NMSSM 


The NMSSM model can solve the p problem in the MSSM and provide a 125 GeV Higgs 
boson without large loop corrections. Throughout this paper we restrict ourselves to the 
NMSSM, where an extended superpotential is given by 



( 1 ) 


After the singlet hied S getting a vacuum expectation value (VEV) an effective /i term 
can be naturally generated 


h — h'e// AUg) 


(2) 
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which is expected to be of 0(10^) GeV. The soft breaking terms in the Higgs sector are given 
by 

Vsoft = XAxSHuHd + ^ ^ 3 ) 

The angle (3 is defined as 

Vu = V sin (3, Vd = V cos (3, (4) 

where Vu and Vd are VEVs of and Hd respectively, and v = \/v'^ + v‘j = 174 GeV. 
Gompared with the MSSM, the tree-level Higgs mass in the NMSSM is enhanced by a new 
term Therefore, a 125 GeV Higgs can be easily achieved in the NMSSM. The 

mass square of the SM-like Higgs is given by 


ml = m| cos^ 2/3 + sin^ 2/3 + 


(5) 


where and denote the loop effects and the mixing effects on the Higgs mass. 

The dynamics of the EWPT are governed by the finite temperature effective Higgs po¬ 
tential, which reads 

Kff = Hpree + Acw + ^CT + Vr -|- Vdaisy (6) 


In Eq. 1^ Vrree denotes the tree-level Higgs potential in the NMSSM. Vcw is the well-known 
Goleman-Weinberg potential at zero temperature H- 


Hew — 




( In 


647r2 




(7) 


where i runs over all particles in the NMSSM, with each having degrees of freedom rij, 
field-dependent mass mi{ipi) and spin Sj. In this work, we adopt the Landau gauge and the 
DR scheme with renormalization scale Q to calculate Hew- In order to maintain tree-level 
relations involving VEVs, counter terms Vct should be introduced. Vr denotes the thermal 
correction at the finite temperature, 

rpi _ pea _ 

Vt = 7 ^ ^ / dx ln(l =F e“ V(8) 

27r Jo 


which can be improved by daisy resummation contributions Vdaisy In practice, this term 
can be achieved by inserting thermal mass contributions in the field-dependent mass. 
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To avoid the baryon asymmetry generated at the EWPT being washed out, the phase 
transition must be strongly Erst order, which can be quantitatively characterized in the 
perturbative calculation as 



( 9 ) 


Here Tc and Vc are the critical temperature and order parameter of the phase transition. 


We implement the effective potential Eq. [^into the public package CosmoTransition |16] 


to evaluate the phase transition numerically and numerically perform a parameter space 
scan. As pointed out in previous works [3^ 021 05], there are mainly three phase transition 
patterns in the NMSSM. When the temperature decreases, the Universe can directly transit 
from the symmetry phase into the electroweak breaking vacuum (Type-III), or undergo 


a intermediate phase in the singlet subspace (Type-I) or the Hu subspace (Type-II). It has 


been found that a strongly phase transition can be achieved in Type-I and Type-III transition 
without light squark contributions to thermal loops |42j . 

The shape of the effective potential at zero temperature, more exactly, the energy gap 
between the symmetry phase and the broken phase (AU = Vsym — Vew) |t=o encodes the 
information on EWPT in the NMSSM. A smaller AV may lead to a lower Tc and thus a 
large transition strength Vc/Tc- This correlation can be understood by 



( 10 ) 


Following this insight, we can perform a semi-analytical analysis of the transition strength 
in terms of the energy gap 02] . 

First of all, the potential energy of the electroweak vacuum can be divided into three 
parts at tree level 


where the individual contributions are given by 


( 11 ) 



( 12 ) 


(13) 


( 14 ) 
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Explicitly, the first part is the contribution from the doublet and is almost hxed by the SM- 
like Higgs mass. The second part is the pure contribution from the singlet. The third part 
results from the doublet-singlet mixing, which can be described by an auxiliary parameter 

c-A, 

C4 = l-(^+^)sm2,3. (15) 

In the case of type-I phase transition, the universe transits into an intermediate phase in 
the singlet subspace before the EWPT. The absolute minimum of the singlet subspace at 
zero temperature locates at the origin or 

Us = (l + a/1 - ) . (16) 

where is given by 

2^, = ^ - g(l + R.f - CA)fvVAl (17) 

where = Ak^/XAi^,. The Tpyp-I transition occurs only if the Ug is the absolute minimum 
with Xg, < 1/9. The relevant energy gap can be easily obtained as |12] 

AH = Vs — Iew ~ ~ 

^ ^ - Ca\^v\uI - vl) + k\vI - ulf 

+ [2ul{us - Vs) + Vs(v^ - uD] , (18) 

where Vs is the tree-level potential energy of Ug. It is obvious to see that a substantial 
deviation of Ug from Vg is crucial to decrease the energy gap away from the doublet limit 
set by the Higgs mass. The numerical results showed that the Type-I phase transition 
correspond two cases, in which either ^ 10 and — 1 < < 0 |12]. For the large 

the deviation of Ug from Vg would be small and is disfavored by the SFOEWPT. 

Note that at tree level the vacua energy gap AH may be negative. Fortunately, the 
Coleman-Weinberg one loop correction would lift Vs and reduce Vew^ thus guarantee that 
the EW vacuum VLew is lower than Vts- However, if the metastable vacuum VLs is well below 
V^EWi the hnite temperature potential could not smooth out a large negative gap when the 
Universe cools down. In this case, the phase transition is not valid. 





Then we consider the Type-Ill phase transition arising in the case that the origin is the 
absolute minimum in the singlet subspace or the origin is metastable. In this scenario, the 
tree-level energy gap is simply given by — Vew as 

+ + (19) 

We can see that the singlet part will dominate the energy gap for a large fi. A small energy 
gap often requires a large negative and —4/3 < < 0. For a moderate p, the mixing 

part becomes important and a negative Ca is useful to decrease the gap. 


III. NMSSM INTERPRETATIONS OF GALACTIC CENTER EXCESS 


The Galactic center excess( GCE) prefers a ~30-40 GeV DM particle annihilating directly 
into bb with a cross-section about (av) ~ 2 x 10“^® cm^/s. In this work, we would like to 
explain this excess by the annihilation of neutralino pair to bb via an s-channel light GP-odd 
pseudo-scalar in the NMSSM. 

The neutralino mass matrix in the NMSSM is written as 


M 


( Ml 


V 


0 

gi'Vd 

V2 

9iVu 

V2 

0 ^ 

M2 

32 

(j2Vu 

V2 

0 


0 

-p 

-\Vu 



0 

-Aud 




2kVs j 


( 20 ) 


The lightest mass eigenstate of the neutralino is the DM candidate, which is dehned as: 


X = 7Vn5 + AisIT + + N,,~S. 


( 21 ) 


The Majorana DM annihilation cross section to bb with a relative velocity Vr for the 
interaction |a(|/axxX7^X + Vathbl^b) is given by 


3 Vaxxygbb 
327rm2 {5 + ’ 


av, 




( 22 ) 
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where 7 = maVal^m^ is a parameter defined by the pseudoscalar decay width Fq, 5 is a 
degeneracy parameter defined as 

2 

777 

If 5 is not very small, the DM annihilation occur off-resonance and is almost velocity inde¬ 
pendent. As discussed in Ref. [ 8 ], the GCE and correct thermal DM relic density can be 
simultaneously accommodated as long as the product yl-^-^Uabb adjusted to an appropriate 
value. In this case, the main components of the DM would be Bino and Higgsino in the limit 
k/X ^ 1. Since the singlet scalar is heavy and is almost decoupled, the SM-like Higgs is 
similar to that in the MSSM. The singlet component of the lightest pseudoscalar can sup¬ 
press its production cross section, hence the stringent limits on the the pseudoscalar mass 
from the LHC H/A ^ searches can be avoided. However, large k/X would enhance 


the — VJvv ^9- (13), which is disfavored by the SFOEWPT. Furthermore, a very small A 
is difficult to give a correct SM-like Higgs mass ~ 125 GeV at tree level. Therefore, in this 
work we do not consider such parameter space. 

In the limit k/A <C 1, the dominant component of the DM may be Singlino. In this case, 
the correlation between k/X and the DM mass in the NMSSM is 

X^v‘^{m^ — /i sin 2/3) 


K 

X 


2/i 




2p’ 


(24) 


m^(m^ — jj?) 

where the second equality corresponds to small The DM would also have non- 

negligible Higgsino components with a moderate fi value, which can be described by for 
small 


iVl3 

iVl5 

Nu 

iVl5 


Xv 

Xv 


cos/3(l 

sin/3(l 








/itan (3 


tan (3) 

)■ 


(25) 

(26) 


For the CP-odd Higgs, the absence of tachyon states favors a negative kA^ or a small 
positive kAk- The constraints from the LHC Higgs researches also compress the allowed 
parameter space. If the heavier CP-odd Higgs is very heavy which is consistent with the 
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Higgs observations, the lightest CP-odd Higgs mass can be written as 


~ ^An^sin2/?(^^^^ -|- 4 k) — 

2 jj, A 


(27) 


If Xv sin 2/3/2/i is not very large, the main component of the lightest psendoscalar is singlet. 
In this case, the conpling gatb is snppressed by a small active psendoscalar fraction. At the 
same time, the conpling //^xx i® determined by small k and small Higgsino fraction. This 
means present DM annihilation cross section shonld be enhanced by the resonance effect 
with small 6. From Eq. we can see that the GCE can be easily explained by adjnsting 
the combination gaxxdabb/i^"^ + 7^) in fhe zero temperatnre limit. In the early Universe, the 
DM annihilation occurring close to the resonance is sensitive to the temperature mm- 
If 5 < 0 and |(5| is not very small, the aVr tends to have a larger value with a larger Vr as 
the temperature increases in the early Universe. In this case, the DM relic density would be 
suppressed. Some non-thermal DM production mechanisms are needed to obtain the correct 
DM relic density. We do not consider such mechanisms and focuss on the thermal freeze-out 


mechanism here. For the 5 > 0 case, from Eq. 22, we can see that the av^ always decreases 
as the temperature increases. This means the process XX a ^ bb can not sufficiently 
reduce the DM abundance in the freeze-out epoch. Some other DM annihilation channels 
are needed to generate the correct DM relic density. 


IV. NUMERICAL RESULTS 

In this section, we perform a numerical analysis to interpret the GCE and SFOEWPT 
in the NMSSM. We use the packages NMSSMTools 4.4.0 [4^ 150] . micrOMEGAs 3.6 .9.2 [5T] 
and CosmoTransition [46] to scan a six-dimensional parameter space (A, k, Ax, A^, /i, tan/3) 
with various experimental limits. Since we focus on the EWPT properties affected by the 
singlet sector, we assume that all the sfermions are heavy and fix soft breaking parameters 
as Ml = Me = 200 GeV and Mq = Mu = Md = At = At = Ai = 2000 GeV. A benchmark 
point, which can simultaneously explain the GGE and SFOEWPT, is given in Tab |TJ 

In order to understand the properties of the parameter distribution favored by the GGE 
and SFOEWPT, we plot Fig. where only two parameters k and vary in the parameter 
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A 

Aa (GeV) 

K, 

A. (GeV) 

tan/3 

p (GeV) 

0.5 

840.0 

0.029 

-99.18 

3.05 

235.0 

m^o (GeV) 

rUa (GeV) 


(cTu) |.^^o(cm3/s) 

(Tsi(cm2) 

PTS 

34.96 

69.54 

0.095 

1.72 X 10“2® 

9.03 X 10“^° 

1.06 


TABLE I: The benchmark point which can simultaneously explain the GCE and SFOEWPT. 
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FIG. 1: Results of a numerical scan with fixed: A = 0.5, Ax = 840 GeV, // = 235 GeV, tan (3 = 3.05, 
Ml = M2j2 = 100 GeV, M^ = 800 GeV, Mj^ = Me = 200 GeV, Mq = M^ = Md = At = 
Afy = Ai = 2000 GeV . The green region denote the parameter space where the DM annihilation 
avr ~ (0.5 — 4) X 10“^® cm^ s“^ can interpret the GCE. In the left panel. The blue region can 
explain the correct DM relic density 0.091 < < 0.138 [521153| via the Z mediated annihilation 

in the early Universe. Also shown are contours of Vc/Tc (the magenta dotted lines) which is usually 
required to be larger than 0.9 by the SFOEWPT. The red and gray regions are excluded by the 
direct detection LUX [55] and the Higss mass 124 GeV< ruh < 128 GeV. In the right panel, we 
show contours of the DM mass (red dotted lines) and mass difference 2m^ — rua (blue dotted 
lines). 


space near the benchmark point and the other parameters are hxed. In the right panel of Fig 
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FIG. 2: Results of a numerical scan with varying A and /r, k = 0.029, = —99.18 GeV and other 

parameters being fixed as Fig. Other discerption are same as Fig. 


red and blue lines denote DM mass and mass difference of 2m^ — rua. The behavior 
of these lines can be easily understood from Eq. 24 and Eq. 27 In Fig we show the 
region consistent with the GCE with an annihilation cross section aVr\vr^o iii fhe range of 


(0.5 ~ 4) X 10“^® cm^ s“^. As expected from Eq. 22, there are two distinct strips near the 
contour 2m^ — rUa = 0 which correspond to the cases 5 > 0 and 5 < 0 respectively. 

If there is no other DM annihilation channel, it is difficult to explain the GCE and correct 
DM relic density simultaneously. In the lower green band for 5 < 0, the DM annihilation 
occurring closer to the resonance in the early Universe would induce a very small relic density, 
while in the upper green band for 5 > 0, the DM particles would be overproduced due to 
the suppression of the DM annihilation cross section in the early Universe. In the left panel 
of Fig the blue region denotes the parameter space where a correct DM relic density can 
be obtained. In this region, the DM annihilation via the exchange of the Z boson plays an 
important role in the early Universe, and sufficiently reduce the overproduction of DM for the 
5 > 0 case. This annihilation channel is a p-wave process and can be negligible in the Galaxy 
with small DM relative velocity. The green and blue bands can overlap at (k, A^) ~ (0.029, 
-99.2 GeV). Since the correct DM relic density requires a non negligible DM-Z interaction. 
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we also check the partial decay width of the Z to a pair of DM which is controlled by the 
Higgsino component of DM. For our benchmark point, this width is 0.59 MeV and is allowed 
by the upper limit ~ 2 MeV from the Z decay measurements [5^ . 

Then we show the parameter distribution of (A, /i) favored by DM results in Fig where 
the other parameters are taken in Tab Similar to Fig there are two bands for the GCE 
correspond to different 6 values. Eqs. and imply that the Higgsino faction of the DM 
and thus the DM coupling to the Z boson are controlled by X/fi. By tuning the value of 
A//i, the correct DM relic density can be obtained via the Z channel annihilation in the blue 
region. The upper green band and blue region overlap at (A,/r) ~ (0.50, 235 GeV). 

In Fig. and we also show the constraint from the direct detection. The LUX 
collaboration has set stringent constraints on the DM with mass ~ (9(10)GeV [33] • The 
spin-independent signature for our benchmark point is dominantly induced by the DM- 
quark scattering process via the exchange of GP-even Higgs. The scattering cross section 
is strongly affected by Higgsino components in the DM. This explains why the parameter 
space with large A and small fi would be excluded by the LUX results. 

We depict the contours of Vc/T^ in Figs and It is can be found that there does 
exist a parameter space satisfying the DM results and SFOEWPT. Eq. implicates the 
relation between the phase transition strength and the energy gap, deeper understanding of 
the SFOEWPT calls for a detail analysis of the energy gap. The energy gap at tree level 
/XVtree and full numerical result with the Goleman-Weinberg one loop correction AVfui are 
shown in Fig. We hnd that the variation tendency of phase transition strength Vc/T^ with 
respect to A(k) and is consistent with those of the energy gaps AVtree and AVfui- 

In our benchmark scenario. Type HI SFOEWPT is preferred since the mixing term with 
a negative Ca tends to sufficiently decrease the energy gap. The parameter dependence of 


the AVtree cau be easily understood by Eq. IV which is rewritten as 

AU ~ ^ sin 2/?HA/in^ + (1 - ^ sin2/3)/i^n^ (28) 


where the hrst term, middle two terms, and last two terms represent the contributions 
from the doublet, mixing, and singlet parts, respectively. For our benchmark model with 
K, ~ 0(0.01) and A,^ ~ —0(100) GeV, a negative contribution of dominates in the 
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singlet part. Hence, AVtree decreases as k and |Hk| increase in the region of k, < —Af^/6vs ~ 
0.035, and then the phase transition is enhanced as shown in Fig. and Fig. [Tj From Eq. 
we also hnd that AVtree is suppressed by the small A in the mixing and singlet parts. The 
dependence of AVtree on is more complicated. In our benchmark scenario with a small k, 
the contribution from the mixing part would be more important than that from the singlet 
part. In this case, it can be found that AVtree becomes more negative as /i decreases when 
p > Ha sin 2/3/4 ~ 110 GeV. 

Note that although the AH at tree level is negative in Fig. the loop corrections from 
the Coleman-Weinberg potential lifts up the AVtree and helps the phase transition to occur 
in the early Universe. The detailed discussions of the loop corrections to AV can be found 
in Ref. [ 12 ]. 



FIG. 3: The energy gaps at tree-level AVtree ( black contours) and the full numerical results 
with loop corrections AVf^i (magenta contours). Left panel: plots of AU in the k-Hk plane with 
parameters set as those in Fig. Right panel: plots of AU in the A-p plane with parameters set 
as those in Fig. 
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V. CONCLUSIONS AND DISCUSSIONS 

In this work, we study the possibility that the GCE and SFOEWPT can be simultaneously 
explained in the NMSSM model. The GCE can be interpreted by the annihilation process 
XX bb via the exchange of a light pseudoscalar. For the Singlino/Higssino DM, the 
annihilation should occur near the resonance with a very small mass difference 2m^—ma > 0. 
In this case, the correct DM relic density is obtained by the DM annihilation via the exchange 
of a Z boson in the early Universe, meanwhile the present DM annihilation cross section 
accounting for the GCE can also be easily achieved. 

Besides the mostly singlet pseudoscalar ~ 0(10) GeV required by the GCE, there is 
also an accompanying light CP-even Higgs in the scalar sector. Such singlet like scalar 
could decrease the vacua energy gap AU with a negative and a moderate p, which 
would lead to SFOEWPT required by successful EWBG. We find that there does exit the 
parameter space in which the GCE and the correct DM relic density can be interpreted by 
the Singlino/Higssino DM with SFOEWPT being realized at the same time. 

The discussions in this work can be extended to the DM models containing new singlet 
states. In such models, the Higgs spectra would be affected by the requirements from the 
EWPT and DM phenomenology. Discovering such new scalar sector at (P (10) GeV is an 
excellent motivation for the future high energy collider experiments. 
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